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Abstract. We study the SU(2) Witten-Reshetikhin-Turaev invariant for the Seifert fibered ho- 
mology spheres with M-exceptional fibers. We show that the WRT invariant can be written in 
\^ • terms of (differential of) the Eichler integrals of modular forms with weight 1/2 and 3/2. By use 

^^ , of nearly modular property of the Eichler integrals we shall obtain asymptotic expansions of the 

^D • WRT invariant in the large-A^ limit. We further reveal that the number of the gauge equivalent 

^^ , classes of flat connections, which dominate the asymptotics of the WRT invariant in A^ — > oo, is 

related to the number of integral lattice points inside the M-dimensional tetrahedron. 
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1. Introduction 

The Witten invariant for the 3-manifold M is defined by the Chem-Simons path integral 
as [58] (see also Ref. 3) 

Zk(M) = fexp(2nik CS(A)) DA (1.1) 

where k eZ, and CS(A) is the Chem-Simons functional 

^, CS(A)=gij J Tr|AAdA + ^AAAAAJ (1.2) 

'nI" ! In a limit fc ^ oo of the Witten invariant Zk{M), we may apply the saddle point method. As the 

\Q ! saddle point of the Chem-Simons functional (11.21) denotes the flat connection 

O 

■> ; dA+AAA = (1.3) 

c^ ■ the asymptotics of the partition function becomes a sum of the Chern-Simons invariants, and it 

is expected to be [3, 1 1, 58] 



^ 



a 

X /j~g-2m(/„/4+dim//0/8)g2m«:+2)CS(A„) q ^x 

Here the sum of a denotes a gauge equivalent class of flat connections, and Tq, and !„ respectively 
denote the Reidemeister torsion and the spectral flow. The first Betti number is Z?^ and H' is the 
cohomology space. 
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To Study the asymptotic behavior of the Witten invariant rigorously, we need explicit expres- 
sion of the invariant. Alternative and combinatorial definition of this quantum invariant was given 
by Reshetikhin and Turaev [46] (see also Ref. 27). We denote Tn(M) as the Witten-Reshetikhin- 
Turaev (WRT) invariant, which is related to the Witten invariant Zk(M) by 

and we have 



Tn 



^ ^ V 2 sin(;r/ 



i(;z-/A^) 

Using this definition of the WRT invariant, asymptotic behavior of the WRT invariants for certain 
3-manifolds has been extensively studied [14,30,31,48-52]. 

Several years ago, Lawrence and Zagier found a connection between the WRT invariant and 
modular form [32]. They showed that the WRT invariant tn{M) for the Poincare homology 
sphere M = S(2, 3, 5) can be regarded as a limiting value of the Eichler integral of vector 
modular form with weight 3/2. Thanks to this correspondence, the exact asymptotic expan- 
sion of the WRT invariant in the large-A'^ limit can be computed, and topological invariants such 
as the Chern-Simons invariant and the Reidemeister torsion can be interpreted from the view- 
point of modular forms. Meanwhile it has been established that this remarkable structure of the 
quantum invariants holds for the WRT invariants for 3-manifolds such as the Brieskom homol- 
ogy spheres [19], 4-exceptional fibered Seifert homology spheres [20], and the spherical Seifert 
manifolds [21]. Also established is a connection between the Eichler integrals of vector modular 
forms with weight 1 /2 and the special values of the colored Jones polynomial for the torus knot 
T,j [24] (see also Refs. 16, 17,25,60) and the torus link7~2,2m [18]. 

One of the benefit of the quantum invariant/modular form correspondence is an observation 
that a limiting value of the Ramanujan mock theta functions [45] in q ^> q'^^^In fj-Qjji outside 
a unit circle coincides with the WRT invariants for the spherical Seifert manifolds [22]. This 
fact opens up a new insight to modular forms and the Ramanujan mock theta functions, and we 
can expect that further studies on the quantum invariant/modular form correspondence should be 
fruitful. 

In this article, as a continuation of Refs. 19,20, we study an exact asymptotic expansion of 
the WRT invariant tn{M) for the M-exceptional fibered Seifert integral homology sphere A\ = 
^{P\,P2, ■ ■ • , Pm), where pj are pairwise coprime positive integers. By use of modular forms 
with half-integral weight, we derive an asymptotic expansion in A^ ^ oo number theoretically. 
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This paper is organized as follows. In Section |2l we review the construction of the WRT in- 
variant for the Seifert fibered homology spheres following Ref . 3 1 . An explicit form of the WRT 
invariant is given. Also discussed is an integral expression of the invariant. In Section|3]we intro- 
duce a family of vector modular forms with half- integral weight. We define the Eichler integrals 
thereof, and study the nearly modular property of a limiting value of the Eichler integrals. By 
use of this quasi modular transformation property, we compute the asymptotic expansion of the 
WRT invariant in the large-A^ limit in Section|4| We shall see that the invariant is a limiting value 
of the holomorphic function [13]. We study a contribution of dominating terms in the large- A'^ 
limit in detail, and reveal a relationship with the number of the integral lattice points inside the 
higher dimensional tetrahedron. Also given is an explicit relationship between the Casson in- 
variant and the first non-trivial coefficient of the Ehrhart polynomial. In Section|5]we give some 
results based on numerical computations. We compare the exact value of the WRT invariant with 
our asymptotic formula. The last section is devoted to conclusion and discussions. 



2. WRT Invariant for Seifert Integral Homology Sphere 

Following Ref. 31, we compute the WRT invariant tn{M) for the Seifert fibered integral ho- 
mology sphere with M-exceptional fibers M = ^(p) = l,(pi, p2, . . . , Pm) where pj are pairwise 
coprime positive integers. Hereafter we use p as M-tuple 

P = (P\,P2,---,Pm) 

The Seifert fibered integral homology sphere l>(p) has a rational surgery description as Fig. [T] 
(see, e.g., Refs. 42, 54, 55), and the fundamental group has a presentation 

h is center 
TTi (XifT)) = L,X2,..., Xm, h I xp = h-^' for 1 < J < M \ (2.1) 

X\ X2' ' ' Xffi = 1 

Here qj e Z is coprime to pj, and we have a constraint so that the fundamental group (12.11) gives 
the homology sphere; 

M 

^y- = i (2.2) 



Here and hereafter we use 



M 
7=1 
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Figure 1 : Surgery description of the Seifert homology sphere I,(p\, . . . , pm) 



When the 3-manifold M is constructed by the rational surgeries Pj/qj on the j-th component 
of n-component link X, it was shown [26, 46] that the SU(2) WRT invariant Tn{M) is given by 



Here the surgery data Pj/qj is encoded by an SL(2; Z) matrix 



;=i 



IjiPj'ij') 



Pj 



The Rademacher O-function 0(t/) is defined by [43] 

' p + s 

Ip r'^ 



O 



q s 



1 



- 12 s(p,q) forqi^O 



for ^ = 



where s(b, a) denotes the Dedekind sum dA.ll) . An nxn matrix L is a linking matrix 

Lj^k = lk( J, k) + — - dj^k 
qj 



(2.4) 



(2.5) 



(2.6) 



where lk(7, k) denotes the linking number of the j- and k-th components of link X, and sign(L) 
denotes a signature of L, i.e., the difference between the number of positive and negative eigen- 
values of L. The polynomial Jki,...,k,X-0 denotes the colored Jone polynomial for link X with 
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color kj for the 7-th component link, and p(U^P'''^) is a representation p of PS L(2; Z) defined by 

yj^Nlql y mod2Nq 

■y=a mod 2N 

for 1 < a,b < N - I [26]. This representation is constructed from 



P(S)a,,= ^J^smi^n 



p(T)a,b = e^" -- 5«,i 



(2.8) 



with 



satisfying 



-ro' -; I 



s^ = (STy = 1 



Based on the fact that the Seifert fibered manifold E(/?) has a surgery description as in Fig. [H 
we have the following result. 

Proposition 1 ( [31]). For the Seifert fibered integral homology sphere with M -exceptional fibers 
M = £(/?), the WRT invariant is given by 



qn[ 4 2J ^e" - Ij Tn{M) 



£i 2PN-1 

e4 



^ , , n«,(e"57"'-e-^") 



Here 4>{p) is defined by 



lyJlPN j;i (e^-i-e-H 



1 ^ IP \ 

<PiP) = 3-- + l2Y^s\^-,Pjj (2.11) 



where s(b, a) denotes the Dedekind sum (lA.ll) . 

Outline of Proof. We use the surgery formula (12.41) . in which the the colored Jones polynomial 
for a variant of Hopf link X depicted in Fig. [T] is given by 

1 n%,sin(^n) 



JkoMu-Mmi-O - 



■&m{nlN) 



h(|;r)] 



M-1 
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Here ^o is a color of component which has a linking number 1 with any other components of £.. 
After some computations using the Gauss sum reciprocity formula (IA.7I) . we get 



e iv I 4 2 



) (e^ - 1) T^{M) 



e-* 



A^-l 



E Z 



1<2PN i^,^. 



a^: _3), 



modpj (e"w'^' — e ~ 



M-2 



*^ 9j (kQ+2Niij)^ I kQ+2Niij . kg+lNiij 

x]~[e"^^^ e^^'"-e"^^'"l (2.12) 



;=i 



We see that the summand in (12.121) is invariant under 

• ^0 ^ ^0 + '2.N and Hj -^ Hj- I for all j, 

• nj -^ nj + pj . 



Using that pj are pairwise coprime, we can then rewrite the multi-sum of (12.121) into a single-sum 
of dTTOb . D 

The WRT invariant can be rewritten in the integral form as follows; 
Proposition 2 ([31]). 

(e^ - 1) T^iM) 



0, N\ A 1) 



Ct 



lyfTPN 



2P-1 



g(z) 



-2ni y Res -^^ + f g{z)dz 

^ z=mN 1 -Q-^'^K J 
C 



m=0 



(2.13) 



where 



g(z) = e-5M- ^ 



.M-2 



and the integration path C passes the origin from (-1 -I- i) 00 to (1 - i) 00. 
Outline of Proof. Key identity is 

2P-1 2PN-\ 

&m{x) = &m{x - N) e"'^'^^ + 2 ;ri y Res huiz, x) + V fuin, x) 

i—i 7=mN i—i 



(2.14) 



(2.15) 



m=0 

Here the function %m{x) is defined by 

®m{x) = I /jm(z, -^) dz 
c 
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where 



hM{z, x) = e" 



.2 



m+2xf;m 



2PN N 



fniz, x) 



I - . - .\M~2 1 _ p-2mz 1 _ fL~2mz 

As the function g(z) is a linear combination of fM(x, z), we obtain the expression (I2.13I) . n 

In view of (I2.13I) . we can decompose the invariant as 

Tf,iM) = T'^\M) + T^^\M) (2.16) 

where t^^^{M) and t^^\M) respectively denote contributions from the residue terms and the in- 
tegral term in (12.131) . It was identified in Ref. 31 that the residue part r)^^(AI) is the contribution 
from irreducible flat connections while the integral term t]^'(A1) is the trivial connection contri- 
bution. The trivial connection contribution is related to the Ohtsuki series [40], and we have as 
follows; 

Proposition 3. In the limit N ^ oo, the trivial connection contribution has the asymptotic ex- 
pansion as 

e-U J(e.-l)r-(M)-);^^^^ (2.17) 

where the T -series is given by 

nl^isinhf— x) 1 ~ T^m 
[sinHPx)]''-^ 2^(2k)l ^ • ^ 

We will discuss later the relationship with the Ohtsuki series. 

Similar integral with T"'(At) in the case of the three exceptional fibers M = 3 appeared in stud- 
ies [1, 36] of the Ramanujan mock theta functions, which still remain to be mysterious and fas- 
cinating topics. This suggests the remarkable fact that the Ramanujan mock theta functions [45] 
are related to the WRT invariant for the Seifert fibered manifolds. See Refs. 21-23 for detail. 

The dominating term of the WRT invariant in a limit N ^> oo follows from the irreducible flat 
connection contributions t^^^{M) as was expected from the saddle point approximation (11.41) . In 
terms of the Witten partition function Zi,(M), the asymptotic expansion of the invariant is given 
by 

M-3 



Zk-2{M) ^ ^ ^M-3"« Z^^\{M) + trivial connection contribution 

domin 
identity 



fl=0 
7(0) 



and the dominating term Z'"^\{M) m N ^ oo can be computed as follows when we use an 



sin(7rz)=7rzf] 1-^ 

n=l ^ ' 
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Proposition 4. In the large N limit, the Witten partition function Zt^{M) for the M -exceptional 
fibered Seifert homology sphere M = l^ip) is dominated by 

2M~2 



_{(B.w\ 2M-3 , 



(M-2)! ^fP 



where Bk{x) is the k-th Bernoulli polynomial (IA.8I) . 



2P-1 r M / 1 

x;^(-ir"s„_,(|i)e-*'- uMji 



(2.19) 



Among a sum of 2 P terms in the right hand side of ( I2.19I ). we can classify the summation by 
the Chem-Simons invariant, which corresponds to an exponential factor -fp mod 1 as will be 
discussed later. 



3. Modular Forms and Eichler Integral 



We introduce the vector modular forms with half-integral weight which play a crucial role in 
analysis of the WRT invariants for the Seifert fibered homology spheres. 

3.1 Vector Modular Forms with Half- Integral Weight 

We set M-tuple 

i ={£,,... Jm) (3.1) 

where £j are integers satisfying < £j < pj. As before we assume that pj are pairwise coprime 
positive integers. For M-tuple £, we define the periodic function ;^^2p(") with modulus 2 P by 

f M 



/ 



XipW = 



( M ^^ 



-We, i{n = P 1 + E^^T 
0, otherwise 



mod 2P 



(3.2) 



where Sj = +1 for Vj. We see that;^2p(") i^ Q^QT^ (resp. odd) when M is even (resp. odd), 

xU-f^) = i-lfxipin) (3.3) 

and that it has a mean value zero, 

2P-1 

J]xipin) = (3.4) 



n=0 



We define an involution aj on M-tuple £ by 

Cj(^) = (^1 > • • • ; ^;-l 5 Pj — ij, ij+\ ,-■ -, im) 



(3.5) 
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for 1 < j < M. As we have 

X2p'^^(f'')=xip(n) (3.6) 

for 1 < i, j < M, the number of the independent periodic functions ;i'2p('^) f^^ given p is 

1 " 
D = Dip) = 2^U (PJ - (3.7) 

7=1 

We note that 

xip{n + P) = -x''J^(n) (3.8) 

for \ < j < M. 
We set 

<? = exp(27riT) (3.9) 



where r is in the upper half plane, r e H. By use of the periodic functions (13.21) we define the 
(^-series by 

^l(r)=^-J]n-^^^^xUn)q^^ (3.10) 



where we mean 



2 

neZ 



l-(-l)^ 

m2(M) = -—— = M mod 2 = 



/', 



when M is even 

(3.11) 

1 when M is odd 



The ^-series 01(t) is proved to be a vector modular form with half-integral weight. The T- 
transformation is trivial, and by use of the Poisson summation formula, 

J] fin) = J] e-2-'"/(0dr (3.12) 

we obtain the transformation formula under the 5 -transformation as follows. 

Proposition 5. The q-series 01(t) is a vector modular form with weight 3/2 (resp. l/2j when M 



is even (resp. odd). Under the S - and T -transformations (12.91 ) we have 

~m2(M) 



<^'^ = \^\ Z^^l^-'/^ (3.13) 



OJ(t+1) = T^O^.(t) (3.14) 



10 
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Here a sum oft' runs over D-dimensional space (13.71) . and matrix elements ofDxD matrices S 
and T are respectively given by 



9M:M-m2(M) P i+v'>^ UlUpY y, LlL 



S' = 



IP 



(-1) 



Y = exp 



fit', M 



n 



sin 



[ Pf 



■n\ 



M 



irv 



t(>-E3^- 



\pi 



(3.15) 



(3.16) 



Our vector modular form 01(t) may be a generalization of modular forms which appear as the 
character of the affine Lie algebra su{l) (a case of M = 1, and ^F^^r) defined below) and as the 
minimal Virasoro model (a case of M = 2) up to the power of the Dedekind rj function. 

For our later use, we introduce other families of vector modular forms. We define the even 
periodic functions by 



e\i{n) = 



1 for n = ±a mod 2 P 
[0 otherwise 

for < a < P, and the odd periodic function by 

±1 forn = +a modlP 
lO otherwise 

for < a < P. We then define two families of ^-series by 

0?V)4Z^>)^^ 



(3.17) 



^g(n) = 



(3.18) 



^f(T) = ^J]n,l,'^^(n)qf 



(3.19) 



(3.20) 



nsZ 



These families are also vector modular forms with half-integral weight. Namely we see that 
0p\r) is a vector modular form with weight 1/2 satisfying 



fo=0 



0^"V+l) = exp| — 7ri|0;V) 



a 



(3.21) 
(3.22) 
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where N is (P + 1) x (P + 1) matrix defined by 

1 



N^ = 



for a = 



IP 



2 lab . ^ 

— cos TT for (3 ^ 0, P 



(3.23) 



1 



y^jip 



cos {b n) for a = P 



j(a) 



Tlie vector modular form *Pp (r) with weight 3/2 fulfills the following transformation formulae; 



v3/2 P-1 



h=l 



¥^"V+l) = exp| — ;ri|¥^"V) 



where Misa(i'- l)x(P- 1) matrix defined by 



/ 2 . lab 



(3.24) 
(3.25) 

(3.26) 



3.2 Eichler Integrals 



The Eichler integral is originally defined for modular forms with integral weight > 2 (see, e.g., 
Ref. 29). In our cases, the vector modular forms Ol(r) have a half-integral weight, so we follow 
a method of Refs. 32, 60 to define a variant of the Eichler integrals. 



We define the Eichler integrals OUr) of the vector modular form Out) by 



O' 



\(j) = Y.^'^ 



m2(M) 



Xipin) q'' 



(3.27) 



n=0 



This can be regarded as a "half-derivative" (resp. "half-integral") of the modular form 01(t) with 

respect to r when M is even (resp. odd). When M is odd, the ^-series 01(t) might be called the 
false theta function a la Rogers [47]. 

Proposition 6. We assume N\ and N2 are coprime integers, and Ni > 0. Limiting values of the 
Eichler integrals <I>1(t) in r -^ N2/N1 are given as follows; 



2PNi 



N2 *2 






k=0 



ipm 



(3.28) 



where Bi,(x) is the k-th Bernoulli polynomial. 
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We note that 
from which we have for N eZ 






N 2P 

= _ (T'^ I nl-m2(M) 

Jt=0 



A proof follows straightforwardly when we use the following lemma [32] (see also Ref. 41). 

Lemma 7. Let Cf(n) is a periodic function with modulus f and mean value zero. Then we have 
as t \0 



i-tr 

i.y-An,^fj 

n=l «=0 

{-ty 

n=\ n=0 

where the L-function is 



^C/n)e-^' - Y,L(-lnXf) ^ (3.29) 

J]nCf(n)e-""' ^YjL(-2n-l,Cf) ^ (3.30) 



^(^'^/) = Z^ (3.31) 



n=l 

/ 



rY,Cf(k)^[s,jj 



Note that the Hurwitz zeta function ^{s,z) defined by 

oo ^ 

^ (n + zr 
has an analytic continuation for k € Z>o as 



(3.32) 



ai-k,z) = -^ (3.33) 

Proposition 8. The limiting values Ou a) of the Eichler integrals with a e Q satisfy a nearly 
modular property. In a limit N ^ oo, we have the transformation formula as an asymptotic 
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expansion as follows; 
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~-z 



L(-2k-l+m2(M),xip) I ni 



k=0 



k\ 



IPN 



(3.34) 



Here N ^X, and a sum of M -tuples £' runs over D-dimensional space. 



Proof We introduce another variant of the Eichler integral by 



^^(^) 



y2n^X (r- 



Oj(T) 



-z (t-z) 



3/2 



dr, when M is even 



%ir) 



1 r'" y^/Tj 
. V2Pi Jz y/r-z 



dr, 



when M is odd 



(3.35) 



where we assume that z is in the lower half-plane, z e H , and z denotes a complex conjugate of 
z. By use of the S -transformation (13.131) . we have 

^ |-ra2(M) 



®:t(^) + 1 T^ 



^sJo^(-l/z) = r^Kz;0) 



(3.36) 



Here r^i{z', a) is an analogue of the period function defined by 

p 



r^ciz; a) = 

p 



Pi 



^iir) 



2n^ X {T-zf 

Oj(T) 



dr, when M is even 



1 r~ 'i'AT) 
——^ — ^^ dr, when M is 

I V2Pi Ja V^-^ 



odd 



(3.37) 



where a e Q. We find that 0^(z) takes the same limiting value with that of 01(t) in a limit 

z, T ^ a 6 Q; 



®^(^) 



= ®^(^) 



(3.38) 



The right hand side of (13.341) arises from an asymptotic expansion of r^i(z; 0), and we ob- 
tain (l334b . ' o 



For our later use, we study differentials of the Eichler integral 01(t), i.e., "fractional deriva- 
tives" of the vector modular form 01(t). From the definition (13.271) of the Eichler integral, we 



14 
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have for b e Z>o 

\ni At) p^ ' ^ 

^ ' n=0 

By the same computation with (13.281) . we have the following. 



2/,+ l-m,(M)^?^(^)^fi 



(3.39) 



Proposition 9. The limiting values of fractional derivative of the vector modular forms ^ It) 
are given by 



(IP d\'~; 



Nj 






N2 ,P- , 



2b + 2-m2(M) ^ 
where Ni > and N2 are coprime integers. 



B 



2h+2~m2{M) 



2PN, 



(3.40) 



The nearly modular property (13.341) of the Eichler integral gives the following asymptotic 
expansion of (13.401) . 

Proposition 10. In the limit N -^ 00, we have the following asymptotic expansion; 



2PN 



N 



■2b+l-m2{M) 



YjXUf^) e^™" B2b+2-m,W (^^) 



n=l 



-1 



12 



-m2(M) 



y i^b+j+i- 



-m2{M) 



j=0 



b-j 



K\ 



(j) 



2Pnj b.moiM) 



X 



2b + 2- m2(M) 
27 + 2-m2(M) 



r 2P 

Zj K"' Zj ^2P(") B2j+2-m2(M) ( J^ 



n=\ 



-f l+E;^UiW 



2b + 2-m2{M)^l^[-'^k-2b-l+m2{M),xip) I ni 



{2P) 



2h+\-m2{M) 



k=0 



k\ 



2PN 



(3.41) 



where the sum of M -tuples £' runs over D-dimensional space, and we have 

b-j-\ 



<=(5)n(i--^- 



(3.42) 
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Proof. We differentiate (I3.36I) with respect to z, and then take a limit z ^ 1 jN. We get 



u — / 



, ;i+m2(M) Zj i' \ dw i^ 

L(-2k-2b-l+m2(M),xip) I n\ 



w^-N 



■ \h oo 



k\ 



IPN 



(3.43) 



In the left hand side, we use 



w 



A \" " At 



dw 



m=\ 



dw" 



\(m) ■ 



where A„ is the Lah number defined by 



m! \m - 1 

This denotes the number of partitions of {1,2,. . .,n} into m lists (a "list" denotes an ordered 
subset) [56,57], and satisfies the recursion relation 

Thus we have 

,h b 



2 dV . 
dw/ 



7=0 



'.m2(M) ^^j 



Here the ^-number is computed as 

which reduces to (13.421) applying (IA.14I) and (IA.17I) . 

As we have (I3.38I) . we get (13.411) with a help of (13. 401) . n 

In the same method, we have formulae of the asymptotic expansions concerning to the periodic 
functions with mean value zero. 

Corollary 11. We assume that C2p(n) is an odd or even periodic function with modulus 2 P, and 
that it satisfies 



• a mean value zero condition, 



C2P(0) = 0. 



2P-1 



2 C2p{n) = 



n=0 
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In the limit N -^ oo, we have the following asymptotic expansions; 
• C2p{n) is odd; 

2PN 

^ ]=Q ^ ' ■' 0=1 c=l 



n=\ 



2 



2b+l ^ L{-2k-2b,C2p) I ^ri 



(2P) 



2b 



Z 

k=0 



k\ 



2PN 



(3.44) 



• C2p(n) is even; 



1 I 1 \ „^A b + 1 



n=l 

^ ;=0 ^ ' ■' «=' 

P 

f=0 

2Z? + 2 'A L(-2A:-2Z?- l,C2p) / ;ri ''' 



2 



{2P) 



2b+l 



k=0 



k\ 



2PN 



(3.45) 



Proof. We use that 



/'p-i 



C2p(n) = ■ 



y C2p{a) <A2p('^) when C2p{n) is odd 

a=l 
P 

y C2p(a) O^'^lin) when C2p(n) is even 



U=i 



In both cases, due to a condition C2p(0) = 0, we see that periodic functions which appear in the 
modular transformation formula such as (13.431) have a mean value zero. Applying Lemma|7l we 
obtain asymptotic expansions. n 
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4. Exact Asymptotic Expansion of the WRT invariants 
4.1 Exact Asymptotic Expansion 

As a preparation to obtain the exact asymptotic expansion of the WRT invariant tn(M) for the 
M-exceptional fibered Seifert homology sphere A1 = ^{p) by use of the vector modular form 
01(t) defined in previous section, we give ^-series identities at the root of unity related to the 
Bernoulli polynomials. 

Lemma 12. We set cofj as the N-th primitive root of unity; 

and assume that a and k are positive integers satisfying < a < N - I. We have 

N-l (a+l)c . ^,^ k r.(y) / / , 1 \\ 

tf(i-<f (^-1)!^ i r'' '[ N II 

where S j^ is the Stirling number of the first kind (IA.14I) . 

Proof. We follow a method in Ref. 2. 
We define the function P(t; k, a) by 

N~\ ^ (a+l)c 

Pit;k,a)=y—^ (4.2) 

where a and k are positive integers. The function P{t; k = 1, a) is computed as follows; 

oo N-l 

Pit;k=l,a) = J]J]aji;'''''-^"'€ 



{a+l)c+nc nt 
n=0 c=l 



1 e^' 



1 - e' 1 - e^^ 



= |(;;rTT^K'('')-^"'«-(i-^))^ («) 



Here we have used 



^ \N-l, if N\n, 



a»«^' — 



£.=1 1-1, Otherwise. 

We further introduce the function Q{t; k, c) by 



Qit;k,c) = - ^ ^, (4.4) 



18 KAZUHIRO HIKAMI 

where k and c are positive integers. By definitions we have 



N-\ 



P(t; k,a) = Y, oji"^'^' Q{t; k, c) (4.5) 

c=l 

The definition (14.41) indicates that the function Q(t, k, c) satisfies a differential-difference equa- 
tion, 

-Q(t- k,c) = k (Qit; k+l,c)- Q(t; k, c)) 
at 

We can check by induction that the function Q{t; k, c) can be written in terms of Q{t; k = 1 , c) as 

*2(''^''^ = (k^. Z(-i)'"^i""' d^2a; he) (4.6) 

m=0 

From (14.51) we find that P{t; k, a) is solved as 

^^^'^'"^ = ^ E(-ir ^r"^ ^A^-^ ha) (4.7) 

m=0 

Substituting (14.31) for the above solution, we complete the proof. n 



Using the arithmetic identity (14.11) . we can rewrite the WRT invariant (12.101) in terms of the 
Bernoulli polynomials. 

Proposition 13. The WRT invariant for the M -exceptional fibered Seifert integral homology 
sphere M = ^{p), which was computed as in (12.101) . is written in terms of the Bernoulli polyno- 
mials as 

1 1 M-3 r.(;+l) 

2m ( iK?) \\ / ln\ \ . 1 1 V~l 1 '^ M 7 

e^(^-.) (e- - 1) r.(M) = -- —— J](-A^)^ -^ 



2PN-1 



/ 1 \ 

^ -i\M 0° «-l N~l M-3 cU) I \2 

^ (-1)^ V V V V V :!iii3 g-f(«-^-^42;^.|) 

CM - 4.V /-/ Aj Zj Zj Zj / 



(M-4)! ^^^^ ^^ ^^^^ 7 

^ ' a=\ b=0 ff(a) n=0 j=l ■> 



^I^'"'^^(^)-^^^<1)| (4.8) 



where for our brevity we have used M -tuple 

(4.9) 
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and y is a signed sum 



M 






We remark that the second term including a sum of a in (14.81) vanishes when 



M . 

T- 



M 

I 

< 1 



iM 1 



Even when y ,_, — > 1, the second term is a finite sum. It is well known that the sum of inverse 
of prime numbers, Yip-, prime ~, diverges, although the sum up to the 10, 000-th prime numbers is 
still 2.709258 • • • . 



Proof of Prof. [721 We first study a case of 2 , — < 1 • In this case, we have 

Pj 

M ^ 2P-1 

-z' Wk^'-^'^'] = Yj^tM)z^ (4.10) 

y'= 1 m=0 

where the periodic function ;^2p("^) is defined in (13.21) . Using this identity in (I2.10I) . we have 
givl 4 2j l^e" - Ij TAr(AI) 

E 2P-1 /V-1 ^n_(,„+(M-3)P)2 2P-1 



Z E^^^« ^^--^ E -"' 



2 VTpTv ;a fe' (eT» - l) 



j=0 



, 2P-1 /V-1 W-1 2m^^ 



m=0 «=0 /t=l ' 



/ 2m* , \'l 

( e A- - II 



. ^ 2P-1 W-l 

M-3 <^0+l) 
X 






Here in the first equality we have used (14.101) . and decomposed a sum of n by setting n = N j + k. 
We have then applied the Gauss sum reciprocity formula (IA.7I) in the second equality, and then 
used our formula (14.11) in the last equality. As we have 

2PA'-l 

Y^/2pin)G^"'=0 (4.11) 



n=0 
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due to (13.41) . we obtain the first term of (I4.8I) . 

For other cases Z; 7" > 1^ the generating function (I4.10I) of the periodic function ;t2p(") ^^ 
replaced with 

M 



X (/(^4-") +(-l)«*i j-'l^'?-)] = ';^'Arf,('>)^" (4.12) 



2P-1 



(1=0 

Thus comparing with a case of 2, — < 1, we need additional term Tadd defined by 

00 fl-l 2PA'-l 21 n- ^i 

£4 e 2PiV 



Tadd - 2j Zj Zj Zj 



^ t^ t^ t^ 2 ^f2PN Un. _ e~^-)^-' 

1^(24-2^-1) ^ (_j^M.l ^-.i^(E4-2.-l)\ ^^_^3^ 



xie"^ 



We decompose a sum of n by setting n = N j + k, and apply the Gauss sum reciprocity for- 
mula (IA.7I) . After some computations, we obtain 



eo «-l N~lN-l Mf(„-^ + iZ,;Hl 

Tadd = 2]^ 2j Zj Zj Zj Zj TT;^ _i,iXM-3 

a=l fo=0 jj(a) *:=! n=() (cv - C w 1 



^ Zj Zj Zj Zj Zj / ._. *\M-3 



Z_J Z_i ' ' ' ' ' ' I 9m 1\ 

a=\ h=0 ifia) n=0 k=l II —Qn] 

where we have used the fact that the sum (14.11) is real. Substituting (14.11) for the above, we find 
that Tadd gives the second term of r.h.s. of (14.81) . and thus we complete the proof. n 



We now aim to relate this expression with limiting values of the Eichler integrals (13.281) and 
differentials (13.401) thereof. For our convention, we introduce an analogue of the Bernoulli poly- 
nomial defined by 

k=m ^ ' 

where M,m eZ satisfying M > m > 0. For a case of m = 0, we set 

/„"W = g^(.+ f) (4.15) 
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Some of explicit forms of the polynomials f^{x) are given below; 






M , 1 

■^' + 240^(5^ + 2) 

Lemma 14. Let the polynomial f^(x) be defined by (I4.14II . Then the polynomial f^-ki^Q^^) '■^ 
even (resp. odd) when k is even (resp. odd). 



Proof. We introduce the generating function of the polynomials f'j^ix) by 

M 

FM(x,y) = Yj'^f^(x)y"'-' (4.16) 

m=0 

Recalling the generating function (IA.14I) of the Stirling number of the first kind, we get 

FM(x,y) = \ [[y + x+--jj (4.17) 



;=i 



which shows that Fm(x, y) is a polynomial of x + y. Furthermore Fm{x, y) becomes an odd (resp. 
even) polynomial oi x+y when M is even (resp. odd). Then we can conclude that the polynomial 
f^_jj(.x) is even (resp. odd) if k is even (resp. odd). n 

By use of the generating function (14.171) . we obtain the following differential equation and 
recursion relation of f^ix); 

^j!:^{x) = {m+\)C{x) (4.18) 



We can rewrite the WRT invariant in Prop.[T3]in terms of these polynomials as follows. 
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Proposition 15. The WRT invariant t^vCAI) for the M-exceptional fibered Seifert homology 
sphere M = 2(p) is written as 



1 1 



M-2 2PN-\ I 



2 (M-3)! 

X e^™ 



2P 



2PN 



+ 



(-ir 

(M-4)! 



oo a-l N-l I „.\l 

y y y'ye-fH-^+i24) 



a=l 6=0 )J(fl) n=0 



-t=0^ I ^ ^ i '^■'J V 



'2n-2Z?-M + 4 + 2/-^ 

J pj 



2N 



l^fri^]BM\ (4.20) 



Proof. To prove for a case of 2y — < 1^ in which the second term in (14.201) vanishes, we only 
need to apply (IA.13I) to 5^+1 (^ ^ ) in the first term of (14.81) . As we have an identity 



J]'-JLBji. + y),J = J]f-(y,-^)B,i.)^ 
j=i J k=0 ^ ^ ' 



(4.21) 



we obtain the required expression. 

In the case of y , — > 1, we need to evaluate the second term in (14.81) . which we have set Tadd 
in the proof of Prop. [131 This term Tadd can be transformed into the above expression (14.201) when 
we apply (14.211) . n 
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Theorem 16. The exact asymptotic expansion of the WRT invariant tj^{M) for the Seifert ho- 
mology sphere M = ^{p) with M -singular fibers in N ^> co is given as follows; 



Ini (<>(?) IW Im 



'I 4 2j ^e" - Ij T/v(A<) 
1 1 



;|-m2(M) Z_i 



M \_l I 1 

2 



2 (M-3)! ii-m2(M) Aj -• \2P;r/ ' l^\,m2(M) 



l^H 



k: 



U) 



X 



1 



2 J + 2 - m2(M) 



r 2P 



n=\ 



-^I l+i:;?7^i 



+ 



1 (-1) 



M 



(M-3)! i 



I -,1 



V V <"'(^), , .M-2/ « m2(M)\ 

Zj Zj^2P («)/2m ^^ —^ 



m=l a=0 



X 



y /^'«+/ 



y=i 



i / 1 

"2 



2P;r 



m-; 



Ky- 



™ 2 - 5c,o - 5c,p „ / c 



K^jrl ^ J] N^ 1Z^5£Z^ B,^{^\ c-^i- 



J 



'2; 



c=() 



2P 



+ 



I M-4 I 

1 f-n^-i L 2 J p 



(M-3)! 

m 



12 



^T^ Y^ <-'(^) M~2 1 a m2(M)\ 



m=0 a=\ 



X 



2/';z- 



/*— 1 

/,-■) 2m + 1 v^ „ / c \ A7c2 

<".I7TtS^°*^'-(27)' ' 

"^ c=l 



1 \ -.0) 2m +1 



+ 



L^J 






i / 1 

2 



ipn 



m-j 



-"-'(PZ, |^-P(2«-i)) 2 - ^,.0 - Sc,P /£_\ t,^, 
J ^ ' 2 'AlPl 



c=Q 

M-4 I 

2 



1 1 Y^^, -1 



+ (-ir^ E ^2'!:;] h 5 - 5 E f L ^ E A^'"--^ 



m=0 



2 2^ Pi i-2 



7=0 



i / 1 

^2 



2P;r 



m-7 



xif' 



(i) 



2m+ 1 



p-i 



m,l 



■' c=l 



'^""^^"--^^""V, (^U--^- 



+ 



E 

/t=0 



r^(fc) / ni 



k\ \2PN 



(4.22) 



w/iere we have used an involution cr on x e Z2P defined by 



cr(x) = P - X mod 2 P 



(4.23) 
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The coefficients Tp{k) in a tail part are defined by 



Tp(k) 



(-ir+i (2^) 



2k 2P 



M-2 



(M-3)! 



Z-r%)Z(-iy^^..(^) 



^_,( n + iM-l)P 
2P 

oo a-l 



n + (M-l)P 



IP 



M-3 



H)0+1) 



] 






xy;' 



M-3 



z> + 



1 1 






( 



B 



2k+j 



1 
JP 



cr 



M-\ 



Py^-(2a-l)P 



(4.24) 



Proof. We first study the case Y^j— <l,'vi\ which we only have the first term in (14.201) . We shift 
the parameter n by P (3 - M), and we have 

xip{n -PiM-3)) = {-If--^-^^^ xT'^\n) 

As we see that;\r?p('^)/f "^(^^^tF^ - [^^^f^J - % is an even (resp. odd) periodic function 
when k is even (resp. odd), we can write by use of f^(x) = -^ that 

2m / ip(p) _ 1 -1 



IV I 4 2, 



(e^ - l) tn{M) 



1 A^- 



M-3 2PiV-l 



2 (M-2)! 



^2^ Xn)e2™"'5M-2 



(_1)M 



1 



n=0 
IPN-l 



2PN 



2 (M-3)! 



Z^- Z^J'%)/ri^ 



+ 



(-1) 



M 



1 



c=l 



r2pw-i 



n m2(M) 



2 (M-3)! 



2P 



eim'^' 52c 



2i'A^ 



c=0 L n=Q ^ ' 



Generating function (14.101) proves 



2P-1 



2^&(")^(") = o 



n=Q 



where g{n) is an arbitrary polynomial of n of order at most M - 1. So we find that 
Xw^^^ fk''^ \2P ~ 2P ~ 2) ^^ ^ periodic function of n with mean value zero, and that the ex- 
pression (14.201) can be identified with a limiting value of the Eichler integrals and their deriva- 
tives studied in the previous section. This proves that the WRT invariant is a limiting value of 
the holomorphic function oi qin q ^> Q2-m/N [\2]. Substituting both (13.441) and (13.451) for above 
expression, we get (14.221) . 
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For the second term Tadd in (14.201) . recalling the periodic functions il/^pin) and Ofpin) we can 
rewrite it into 



s^M °° Q~l 



'^'"^'^ ~ 2(M-4)! ) ' ) ^ ) ' "! ) ^ ^^^ ^ ^2m ^ 



L^J 



2 1^11: 



Pi) 



zzz z^-/^^ 

a=\ b=Q rf{a) ^ m=0 

2PN-1 

Z(o.M-l(P2,f-(2a-l)P)) . „2 I n \ 



n=0 



L^J 



+i-ir-^j]N-"fZi]\b+'--'-j]l 



m=Q V 

X 



/ ^'J 



/i=0 



((r'^-'(P2,^-(2a-l)P)) 



2P 



(n)e'^^^ 52^+1 



2/'A^ 



With this term, the mean value zero condition is satisfied even in this case, and we can apply the 
result of Coro. [TT]to obtain the exact asymptotic expansion as required. n 



We have thus obtained that the WRT invariant Z/v_2(A1) for the M-exceptional fibered Seifert 
homology sphere M = ^{p) is written in a limit A'^ ^ oo as a sum of exponentially divergent 
terms and a tail part; 



M-3 



Zn^2(M) ^ Yj ^"'^''^ 4^-2(^) + a tail part 



(4.25) 



k=0 



Here a tail part means an infinite power series of I IN, and it corresponds to a contribution from 
the trivial connection. Among the divergent terms in (I4.25I) . the dominating term in the limit 

N ^ oo\^ Z''^\^{M), which is read as follows. 

Corollary 17. In a limit N ^ oo, the asymptotics of the WRT invariant TN(Ai) for the M- 
exceptional fibered Seifert homology sphere M = ^{p) is dominated by N'^~^ ■ Z^_2(A1), namely 



Zn-2{M) ~ A^^-' • e"^ 



•HP) ■ 1 



;m2(M)-l p-4m 



e 4" 



2 V2(M-2)! 



X 



2sf'^^^C,(^e-(^^^-^^". (4.26) 
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I.e., 



)M-2 



Z^_2(M) ~ A^' 



M-3 



(M-2)! yfP 






Xc,(f>-('<-^) 



-i- mW 



^ -,^Mp('l+v,■■^Wpv,■-!-+py,■Vw,■-^ 

X (— 1) I ' I'il ^ I'j '^1 '^''*i I'jPk 



M 

n 



sin 



pf J 



(4.27) 



where the sum of M-tuples t runs over D -dimensional space (13.71) . and the function Cp (ij is 
defined by 

IP 
Cp [^ = Yj^ipin) Bm-2 (^) (4.28) 



n=\ 



We note that the invariance (13.61) of the periodic function ;^f2p(") indicates that 

Cp(aiaj{i)) = Cp[i) 



(4.29) 



By construction, the asymptotics (14.261) should coincide with ( 12.191) which follows from 
residue part of (12.161) . We do not have a direct proof, and we have checked the equivalence 
numerically for several jtfs. Recalling the path integral approach, the sum of M-tuple /can be 
regarded as a label of the gauge equivalent class of flat connections a in (11.41) . and we can identify 
the Chern-Simons invariant with 

.2 



M o \ 



ij 



csK,) = -^ i.E^ 



-iPn 



mod 1 



(4.30) 



See Refs. 10, 28 for computations of the Chern-Simons invariant for the Seifert homology 
spheres. Note that this value originates from the T-matrix (13.161) of the vector modular form. 
Correspondingly, the Reidemeister torsion is given by 



Ta(t) - 



M / 

sm 

j=i 



n 






Cp(^ 



(4.31) 



Here the product of sin-functions originates from the S -matrix of the vector modular form. 



4.2 Ohtsuki Series 



A tail part in the asymptotic formula (14.221) has a simple generating function as was studied 
in dTTTt . 
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Theorem 18. Let the T -series be defined by (14 .241) . Then the generating function of the T -series 



IS 






Proof. We first study a case of Xiy 7- < 1- Using (I4.10I) . we have 



SSP^^II'f^M^^— "«' 



We equate this expression with ^^q jj^ x^^. Applying the Mellin transformation, we get 

M-3 



(_1)M (2Py'' Y^' g "<-< 1 (n + P(M-3) 

'^'-~r W^. ^ ^^"^"^ 4^ TTJTik ^'^^''' 2P — 

^ ' n=0 j=Q ■' ^ 



— ik) (k\(PiM-3)-n''-^ 



xZ^M-3 



, 2P 

k=J 



Here we have used (13.331) . Identities (14.141) and (14.191 ) give 



V.(.) '^^ '/'(M-3)-.^^^ l)--'(/il)r-^(" '^ 


en we have 




^ 1 (2P)2^ ^^^ . 




X 




3r this expression, we substitute an identity 


(n + P{M-3)\ 


I n m2(M) - 1 \ 
^2/' ' 2 j 




+ {2k+ ) y -—+m+ ' 
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which follows from (lA.lOl) . As we have 



M-2 



Zifr 



y=i 



n \\ I n m2(M) - 1 

IP 2/1 IP 2 



7-1 



r, i n I n m2(M) - 1 ' 
= pM-2\-^-^,-j^-m ^ 1 = 



from (14.171) . we obtain (I4.32I) . 



When Zj ~ > 1, we have another term coming from the second term in (I4.12I) . This gives an 
additional term to (14.331) : 



'^add 



(-ir (2^) 



2k oo a-\ 



2 (M - 4) 



oo a-\ M-3 . f -. -. ^ 



fl=l b=0 if{a) j=l 



2k + j 



x\B 



2k+j 



M-A 1 



2 Zj p. 



+ (-iy5: 



ik+j 



M-2 1 



^ ,■ Pi J) 



Recalling Lemma[T4land applying the same method with above, we recover (14.241) . 



D 



We give explicit forms of some T-series as follows; 
Tp<0) = 
Tpi\) = 4P 



T,(2) = SP' 



TA3) = AP' 






5 Zt^+^-M -2 Z7T + 2-M 



,i=i Pi 



W=i P.- 



// 



(4.34) 
(4.35) 

(4.36) 
(4.37) 



Based on the exact asymptotic expansion (14.221) of the WRT invariant, we extract the tail part 
and define the formal ^-series roo(AI) as a quantum invariant of the Seifert homology sphere 
M. = 'L(p) by identifying exp(2ni/N) with q; 

(4.38) 



M_i , ,, ,.., yr^Tp{k) l\ogq\ 

q^ K^ - 1) • T«,(M) = ^ ^p 1^ 



This invariant Too{M) of the formal ^-series coincides with the invariant T](J'(At) defined in (12.161) . 
Namely we have an integral expression for Too(AI). 
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The Ohtsuki series [40] An{M) is defined from the formal ^-series Too(AI) by 

oo 

T^{M) = J]A„(M)(q-ir 

n=0 

Then the Ohtsuki series An(M) for M = ^{p) is computed as follows. 
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(4.39) 



Theorem 19. The Ohtsuki series A„{M)for the Seifert homology sphere M = ^(p) is written in 
terms ofp and (pip) defined in (12.111) as 



An(M) 



where the function G{x) is 



(n+ 1)! 



^ /p d^ 1 (f>(^ ; 



L7=0 



1 luHr2 ^ 9 



G(.x) = 



n;isinh(^) 

[sinh(;c)]'^"^ 



-J 



G(x) 



x=0 



(4.40) 



(4.41) 



Proof. From (14.381) and (IA.15I) . we have 



m=0 ^ ' y=0 yt=0 •' ■ ^ ' 



We then have 



^ 1 (L-mx^ 



;=o 



(j+1)! \ n-j 



k=\ 



(k) Tp{k) 



n n+ 1 -m 



~ (n+ 1)! ^ 2-i 



m + fc\ .(,„,,) /I 0(/?)r r#) 



(«+l)!j^o ^ \ "^ /""' \2 4 / (4P)^ 



where in the second equality we have expanded the binomial coefficient in terms of the Stirling 
number of the first kind using (IA.14I) . then we have applied (IA.17I) . Theorem fT9l shows that the 
function G{x) (14.411) gives 



, 2/t 






x=0 



i^Tpik) 



Substituting this expression and recalling (IA.14I) . we obtain the required formula. 
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Explicit forms of the lowest 3 Ohtsuki series An(M) for M = l^{p) are 
UM) = 1 



{AM) 
(4.43) 



A2{M) 



3 {(f>(p)f + 12 cl>(if)- 4 P ( 



96 



16 



M 






i^iP) + 2) 



p2 

"96 






where AdM) is the Casson invariant of the Seifert homology sphere M = l^{p) [12, 39] 



Ac{M) = -- + 



8 24P 



hZ(£f-<--«^-5§^(£''") 



(4.44) 



(4.45) 



k=\ '^' '^' J k=i 

The relationship between the Casson invariant AdM) and the Ohtsuki series Ai(M) was first 
proved in Ref. 38. See Ref. 53 for a computation of A2iM). 



4.3 Lattice Points 



We have seen that the asymptotics of the WRT invariant is dominated by the term (14.261) . 
which shows that the number of terms in the sum of £ are at most D defined in (13.71) . Though, 
as was studied in Refs. 19,20 for cases of Af = 3 and M = 4, the function Cp(i) may vanish for 
some ^'s. 

Theorem 20. We fix M -tuple p with pairwise coprime positive integers pj, and let the function 
Cpit) be defined by (RTIHl for i ^ 'L" satisfying 1 < £j < pj - I. Due to K29[i . we have D 
independent functions. We set y{p) as the number of M -tuples t satisfying 

Cp{€) ^ 0, 

and L(p) as the integral lattice points £ e Z^ inside the M -dimensional tetrahedron, 



Then we have 



o<y^<i 



D-y{p)> L{f) 



(4.46) 



Proof. As a generalization of the function Cp{i) defined in (14. 281) . we define CKii) for fc > by 



2P 



c'p{e) = Y,xUn)BA-^ 



n=\ 



(4.47) 
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in terms of the Bernoulli polynomials. We have Cp(£) = C^~^{£). As a generating function Z^(?) 
of these polynomials, we define 

k=0 



Using (IA.8I) . we have 



4^^) = iTTT Z^2p(«) e*" (4.49) 

n=l 

In the case of < Z , — < 1 , we have 

^y Pj ' 

n=0 ;=1 ^ ' 



which gives 






This shows that 



zi{t) = - 






^M ^ ^(-^M+2^ 



and that Cl(^) = for < yt < M - 1 . So we have Cpi£) = when < 2; ^ < 1- 

The invariance (13.61) of the periodic functions ^jpC") proves the statement of the theorem, n 

In the case of 2/ — > 1^ the generating function (14.501) should be replaced with a formula 
like (I4.12I) . and we do not know whether the function Cp(£) vanishes. We conjecture that, when 
2^- — > 1, we have Cp{£) = iff^jpC") coincides with;^f2'p(") S-t. Zy 7" < 1- 

Conjecture 1. Under the conditions of Theorem\20\ we have 

D-y(p) = Lip) (4.51) 

This conjecture was proved for M < 4 in Refs. 19, 20. It states that the number of the flat 
connections which contribute as (14.261) coincides with the number of integral lattice points inside 
the M-dimensional tetrahedron. 
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4.4 Ehrhart Polynomial 

Explicit form of the number L(p) of the lattice points inside the M-dimensional tetrahedron 
was first computed by Mordell [37] for cases of M = 3 and M = 4; 

• M = 3; 

L(pi,P2,P3) = ^(Pl- 1) (P2 - 1) (P3 - 1) + J2p ~ 4 

-77:1 2 2 2~ ■^(Pl ^2, Py) - S(P2 P3, Pi) - S(P1 P3,P2) (4.52) 

12 \ p^^ p^ p^'j 



• M = 4, 

L{PuP2,P3,Pa) 

1 rr 3 P P v% 1 + p, 1 'f^ 



4 or.r,4i,„ -,/' 4^ 

8 1 r^^ '^ ■ 8 12 ■ 24 ^ p.2 ■ 24/^ ' ^^^ 

4 , 1 4 / „ \ , 4 



-45^4|:.(j..|45.(^..) »... 



Here ^(Z?, a) is the Dedekind sum (lA.ll) . For higher dimension M, the lattice points L(7?) might 
be written in terms of Zagier's higher-dimensional Dedekind sum [59], but there seems to exist 
no applicable expressions. 

Although, there is a useful tool to count the lattice points (see, e.g., Ref. 6). 
Let P be the M-dimensional open tetrahedron with integer vertices, (pi,0, . . . ,0), 
(0, 7?2, 0, . . . , 0), . . . , (0, . . . , 0, pm\ and (0, . . . , 0); 

{ \ " i 1 

^= (^i,...,^M)eZ^ y^<l,4>0 (4.54) 

Let Ep{t) denote the number of lattice points in the dilated tetrahedron tV. So we have 

L{p) = Ep{t = 1) 
In the same manner, we suppose that Ep{t) denotes the number of lattice points of the closure of 

£^0 = #|(mi,...,mM)6Z^| Yj'^<t,mk>o\ 

These functions, E-p(t) and Ep{t), become polynomials of t [9], which are called the Ehrhart 
polynomial. Moreover we have the Ehrhart-Macdonald reciprocity formula [9, 33], 

E^(-t) = (-1)^ E^(t) (4.55) 
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In general, the number of lattice points Ep(t) becomes polynomial of t for arbitrary polytope 
f [9]. We set the coefficients of the Ehrhart polynomial as 

Ep{t) = Cm{P) t"" + Cm-i(^) t""-' + • • • + co(^) (4.56) 

It is well known that cm(!P) is the volume off, cm(^) = Vol(f), CAf-i(f )is ahalf of the boundary 
surface area, Cm-\{P) = j Vol(5f). The coefficient Co(P) is the Euler characteristic ;^^(!P), and 
cqCP) = 1 when P is the convex polytope. 

The first nontrivial coefficient of the Ehrhart polynomial for the M-dimensional tetrahedron is 
thus CM'iiP)- In our case of the M-dimensional tetrahedron with pairwise coprime integers pj, 
we have [5, 8] 



(M-iy.-CM-iCP) 



_ M 1 



/ M / „ \2 r M ^\^\ 



4 24P 



k=i ^P^ 



P 



P 



\k=i P^J 



M , p 



El^h^E^ -Z^l^-) <^-'^' 



k=\ 



Recalling the Casson invariant Ac{A\) (14.451) for the Seifert homology sphere M = 'L{p), which 
is proportional to the first Ohtsuki series (14.431) . we have the following. 

Proposition 21. The Casson invariant Ac{M) for M = ^(p) is related to the first nontrivial 
coefficient of the Ehrhart polynomial for the M-dimensional tetrahedron P (I4.54I) ." 



AdM) - ^ ' cm-2{P) = S— + -^^^-o X (4.58) 

7 X 74 X ^—i n.- rti. 

\<j<k<M 



It is remarked that the residue formula for Cm~2{P) given in Ref. 5 looks like an expres- 
sion (ITTOl) of the WRT invariant tn{M) for M = S(^. 

5. Some Examples of Numerical Experiments 

We give some numerical experiments on the asymptotic behavior of the WRT invariants for 
the Seifert homology spheres. 

5.1 2(2,3,5,7,11) 

For p = (2,3,5,7, 11), we have P = 2310, D = 30, and (pip) = ^. The bases for 30- 
dimensional space is given by ^ = (1, 1,^3,^4,^5) with 1 < /"a < 2, 1 < ^4 < 3, 1 < /"s < 5. 
For all these 5-tuples £, we can check that Cp(£) + 0, which supports Conjecture [T] as we have 

y.± = 2927 ;^ 
^1 Pj 2310 • 
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A^ 


exact result for Z^ 


asymptotics Z^^ 


22 


-13.346013 


+ 17.3979061 


-12.2403 + 16.70131 


23 


-0.57682556 - 


0.511081471 


0.020572 + 0.0041401 


98 


0.93263590 - 


0.49655457 i 


0.323366 + 0.00570231 


99 


22.826764 


-367.893601 


22.8460-365.8701 


100 


464.33437 


-287.595561 


475.688-287.9731 


998 


9.2292110 


-9.33241291 


10.7013-1.605811 


999 


-52995.123 


- 87204.076 i 


-53072.7-87187.81 


1000 


694.74344 


+ 9181.29351 


683.369 + 9183.491 


2398 


-64.891808 


+ 46.6207941 


-62.4971+47.92751 


2399 


320910.08 


+ 27551.3951 


321128. + 27510.11 


2400 


142206.21 


-1871.80801 


142145. -1869.061 


2401 


214250.48 


-80025.1871 


214270. -79907.41 



Table 1 : The WRT invariant Z^iM) for M = 2(2, 3, 5, 7, 1 1). Asymptotic formula for Z^^\M) is from 622 



In tabled we give numerical results on the Wltten Invariant Za?(A1) for M. = 2(2, 3,5,7,1 1), 
which Is performed with a help of PARI/GP. We give both the exact value Za?( A1) and asymptotic 
value Z^^AI) for several A^^s. They vary much with A^^, and comparing these data we see an 
agreement. 



338099 
408408 



and D = 5040. This D = 5040-dlmenslonal 



5.2 1(3,7,8,11,13,17) 

For ^=(3,7,8,11,13, 17), we have (p(p) 
vector space Is spanned by / = (1, 1 < ^ < 3, 1 < 4 < 7, 1 < ^ < 5, 1 < ^5 < 6, 1 < 4 < 8). 
Among these D = 5040 bases, we check that C^(/) = when /= (1, 1, 1, 1, 1, 1), (1, 1, 1, 1, 1,2), 
(1,1,1,1,1,3), (1,1,1,1,2,1), (1,1,1,1,2,2), (1,1,1,1,3,1), (1,1,1,2,1,1), (1,1,1,2,1,2), 
(1, 1, 1,2,2, 1), (1, 1,2, 1, 1, 1), (1,2, 1, 1, 1, 1), which supports Conjecture[T] 

Numerical results on the exact value and the asymptotics of the Wltten invariant for Z(;?) are 
summarized In Table |2| We see an agreement. 

6. Conclusion and Discussion 



We have studied the asymptotic expansion of the SU(2) WRT Invariant tn{M) for the M- 
exceptlonal fibered Seifert homology spheres M = ^{p) in A/^ ^ oo number theoretically. We 
have found that the invariant can be written in terms of a limiting value of fractional derivative, 
i. e. derivative of the Eichler integral, of the vector modular forms with weight 3/2 and 1 /2. This 
supports a result [13] that the WRT invariant is a limiting value of the holomorphic function in a 
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A^ 



58 

59 

60 

61 

118 

119 

120 

121 

238 

239 

240 

241 

242 

243 

244 

998 

999 

1000 



exact value Zf^{A\) 



365.32895 + 679.07006 

1331.8460-433.95047 

-944.99493+765.34451 

130.91099 + 2814.5744 

-0.8206017 + 61.590246 

8.1857781 + 13.369868 

5259.2853 + 4064.4029 

8733.0140 + 5274.8273 

-219.36738-1.608943 

-6151.0562-5617.75586 

-11.492746 + 6.1192358 

-26057.019-52201.108 

49736.853 - 46390.033 

189895.62 + 265408.04 

3782.8814-12474.142 

21039.448 + 18091.568 

-12.505553+49.861847 

78229.306 - 164203.36 



asymptotics Z^-1\M) 



351.149 + 691.982 
1358.51-437.953 

-915.949 + 742.606 
62.8489 + 2763.93 

0.782372 + 60.1248 

0.0195662 + 0.0062675 

5232.38 + 4043.94 

8659.21+5338.15 

-216.499 + 1.53462 

-6220.64 - 5620.95 
1.67454 + 2.34920 

-25950.5 - 52634.8 
49818.0-46337.0 
189029. + 265225. 
3814.35-12433.5 
21107.1 + 18191.2 
-0.0331549 + 0.0338852 
78333.1-164618. 



Table 2: The WRT invariant Zn{M) for M = 2(3, 7, 8, 1 1, 13, 17). Asymptotic formula for Z^^\M) is from (g^l 



limit that q tends to the A'^-th root of unity. By use of the nearly modular property of the Eichler 
integral, we have obtained an asymptotic expansion of tn{M) in the large A^ limit. 

Although an asymptotic behavior of the WRT invariant was previously studied in, e.g., Ref. 31, 
the correspondence between modular forms and the quantum invariants enables us to relate topo- 
logical invariants such as the Chern-Simons invariants, the Reidemeister torsion, and the Casson 
invariant, with geometries of modular forms. For example, we have found that the number of the 
gauge equivalent classes of flat connections, which dominate the WRT invariant in the large-A/^ 
limit, is related to the number of integral lattice points inside the M-dimensional tetrahedron. 
From this view, we have established that the Casson invariant for the Seifert homology sphere 
has a relationship with the first non-trivial coefficient of the Ehrhart polynomial. 

In our previous papers [21, 22] we have shown that the WRT invariants for the Seifert man- 
ifolds with 3-exceptional fibers coincides with a limiting value of the Ramanujan mock theta 
functions. Investigated [23] is a modular transformation formula of the newly proposed mock 
theta functions based on explicit form of the WRT invariants. This intriguing correspondence 
seems to originate from a result that the integral expression (12.161) of the WRT invariant has a 
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connection with the Mordell integral [36]. Our results presented here will shed a new light on 
geometric and topological aspects of modular forms. 

Though we have studied only the SU(2) invariant, the Witten partition function (11.11) can be 
defined for arbitrary gauge group, and an explicit form of the invariant for the Seifert manifold 
is given [15]. Extending the method of Lawrence and Rozansky, it is shown [35] that the parti- 
tion function can be written in the integral form which can be interpreted as the matrix model, 
and that it becomes a sum of local contributions from the flat connections. This fact is recently 
reinterpreted from the viewpoint of the path integral by use of non-abelian localization [4]. As 
it is well known that the Chem-Simons perturbation theory of the SU(AO Witten invariant as a 
\/N expansion can be interpreted from the string theory (see, e.g., Ref. 34), it will be interest- 
ing to investigate the quantum invariants/modular forms correspondence for the WRT invariant 
associated with SU(A^) gauge group as a generalization of the present work. 
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A. Special Functions and Identities 
A.l Dedekind sum 

The Dedekind sum is defined by (see, e.g., Ref. 44) 

*'-) = E((5))((")) (A.1) 



yt=l 

where ((x)) is the sawtooth function 



m) = 



x-lx} - - when x i 



when a: e Z 

The Dedekind sum can also be written as 



1 ^ Ik \ Ikb \ 
s(b,a) = -— y cot -;?r cot — n] (A.l) 

4«^ \a I \a I 



The Dedekind sum is known to satisfy the reciprocity formula 



1 I la b 1 , 

s{b, a) + s{a,b) = -- + —[- + - + —■] (A.3) 

4 12 \o a ab 
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We note 

s{-b, a) = -s{b, a) (A.4) 

s(b, a) = s(c, a) for be = I mod a (A. 5) 

A.2 Gauss sum 

As a discrete analogue of the Gaussian integral, we have a formula of the Gauss sum as 

2^-1 
n=0 

The reciprocity formula of the Gauss sum follows from the Gauss integral as (see e.g. Refs. 7,26) 

y g^f n2+2;rito ^ / — e^ sign(WM) y ^-ni^in+kf ^j^j^^ 

n mod N n mod M 

where N,M eZ and Nk eZ, and A^ M is even. 

A.3 Bernoulli Polynomial 

The n-th Bernoulli polynomial B^ix) is defined from the generating function as 

;7TT = Z^^«^ (A.8) 

k=Q 

Some of them are written as follows; 

Bo(x) = 1 

1 
Bi{x) = X- 2 

B2ix) = x^ - X + - 
6 

,3,1 
B3(x) = x' --x^ + -X 

These polynomials satisfy the following relations; 

B,(\ - X) = (-\f B,(x) (A.9) 

Bkix + 1) - Bkix) = kx^-^ (A. 10) 

—Bn(x) = nBn^,(x) (A. 11) 

ax 
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Note that the Bernoulli function has the Fourier expansion as 



and that 



Bk(x - UJ) = -kl y —-—-^ (A.12) 

Bnix + y) = J]hB,(x)y'''' (A.13) 

k=o ^ ' 



A.4 Stirling Number 



The Stirling number of the first kind 5jf' ^ denotes the (signed) number of permutations of n 
elements which contain m permutation cycles (see, e.g., Ref. 57). The generating function of 
Sj,™^ is written as 



n-l 



f[(x-j) = ^5rV (A. 14) 



7=0 m=0 

. (m) 



and Sn i" when n > m>0. Another form of the generating function is given by 



m! 



E^r^^ (A.15) 



(m) 



Based on these generating functions, we have the recursion relations of 5„ as follows; 



^r= E(j^ii^r" (A.17) 
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